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Fractional topological phases in three-dimensional coupled-wire systems 
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It is shown that three-dimensional systems of coupled quantum wires support fractional topo¬ 
logical phases composed of closed loops and open planes of two-dimensional fractional quantum 
Hall subsystems. These phases have topologically protected edge states, and are separated by ex¬ 
otic quantum phase transitions corresponding to a rearrangement of fractional quantum Hall edge 
modes. Some support for the existence of an extended exotic critical phase separating the bulk 
gapped fractional topological phases is given. Without electron-electron interactions, similar but 
unfractionalized bulk gapped phases based on coupled integer quantum Hall states exist. They are 
separated by an extended critical Weyl semimetal phase. 
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I. INTRODUCTION 

Since the experimental discovery and theoretical ex¬ 
planation of the fractional quantum Hall effect frac- 
tionalization in interacting topological systems has been 
an important theme in condensed matter physics. While 
the understanding of fractionalized topological phases 
has impressively developed in two dimensions (2D), 
much less is known in three dimensions (3D). Slave- 
particle approaches have for instance allowed one to make 
progress for topological Mott insulatorsand 3D frac¬ 
tional topological insulators A”— Some exactly solvable 
models have been reported /^d^ and the Kitaev honey¬ 
comb modeli^ has been generalized to 3D lattices 
Further examples for fractionalized 3D phases include 
spin ice^ — and stacks of 2D fractional quantum Hall 
layersIn the latter, inter layer couplings can stabi¬ 
lize many-layer versions of Halperin bilayer states^”— 
or exotic phases with fractionally charged fermionic 3D 
quasiparticles Other coupled-layer construction s^^d^ 
can also exhibit string-like excitations, which in general 
allow for non-trivial 3D braiding 

Instead of coupling extended layers, this work engi¬ 
neers 3D fractional topological states by connecting 2D 
building blocks of finite width along different directions. 
While coupled-layer physics can be recovered by dom¬ 
inantly coupling the building blocks along planes, the 
blocks may also connect along other geometries, such as 
closed loops. This gives rise to additional topological 
phases and phase transitions. For concreteness, the re¬ 
mainder studies narrow integer and fractional quantum 
Hall strips as building blocks. As a further difference to 
previous studies, the individual building blocks are con¬ 
structed from coupled quantum wires containing spin- 
polarized electrons. Time-reversal symmetry is thus bro¬ 
ken from the outset. Since each wire can be treated as 
a Luttinger liquid^ this approach is especially powerful 
for the analysis of interacting phases. Starting with the 
pioneering work of Kane et it has been shown 

that coupled-wire constructions allow for an analytically 
tractable description of integer and fractional. Abelian 
and non-Abelian topological 2D statesCoupled- 


wire constructions have also led to qualitatively new re¬ 
sults, including a classification of interacting topologi¬ 
cal phases^ and the prediction of spontaneously time- 
reversal-symmetry-broken states towards which 2D frac¬ 
tional topological insulators can be unstable^^ 

The present work extends this list by adapting coupled- 
wire constructions as a tool for the analysis of novel in¬ 
teracting topological physics in 3D. The potential of this 
approach is exemplified by constructing a system that 
hosts several fractional topological phases. The coupled- 
wire construction provides simple illustrations of these 
phases in terms of closed loops (or cylinders) and open 
planes of integer and fractional quantum Hall subsys¬ 
tems. It also allows one to identify a regime in which 
an extended exotic critical phase may exist. At special 
Luther-Emery-type points, the low-energy theory of the 
system can furthermore be solved by an exact mapping 
to non-interacting fermions of fractional charge. 

The plan of the paper is as follows. Section UTl discusses 
the individual building blocks that are used to construct 
the 3D system. Section [111] is devoted to the analysis 
of the non-interacting array, whose phase diagram is de¬ 
tailed in Sec. lYl In Sec. 01 turn to the description of 
the interacting 3D system in terms of coupled Luttinger 
liquids. The interacting phase diagram is discussed in 
Sec. [Vll The results are finally summarized in Sec. IVIIl 

II. INTEGER AND FRACTIONAL QUANTUM 
HALL BUILDING BLOCKS 

The Abelian, but in general fractionalized, phases 
studied in the remainder are constructed from two fun¬ 
damental building blocks A-X-B and C-Y-D shown in 
Fig.E The dispersions of the electron-type wires A and 
B, and the hole-type wires C and D, are asymmetrical 
with respect to zero momentum along the wires. For 
a given chemical potential /i, the Fermi points reside at 
momentum kx = —ki ^2 in wires A and C, and kx = -\-ki ^2 
in wires B and D. This can be realized in spinful wires 
with spin-orbit coupling, which are polarized by a mag¬ 
netic field parallel to the spin-orbit direction. The cen¬ 
tral wire X is of hole-type, while Y is of electron-type. 
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FIG. 1: The building blocks A-X-B and C-Y-D, and the dis¬ 
persions E^'\kx) of the different wires. On the left, strong 
(weak) tunneling couplings are indicated by thick (thin) lines. 


They have Fermi points at kx = In each build¬ 

ing block, the close-by inner and outer wires are tunnel- 
coupled by a strong hopping term. A direct tunneling, 
albeit of reduced strength, also exists between the more 
distant outer wires. 
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FIG. 2: Section of the periodic 3D system (the dotted box 
shows the unit cell). Thick dotted diagonal lines depict the 
dominant hoppings in each A-X-B and G-Y-D building block, 
thin dotted diagonal lines indicate subleading hoppings be¬ 
tween the X and G, D (Y and A, B) wires. The edge mode 
couplings are tyi along solid lines, ty 2 along dashed lines, tzi 
along dotted vertical lines, and tz 2 along dash-dotted lines. 


III. 3D SYSTEM OF NON-INTERACTING 
WIRES 


Reference m showed that the building blocks form nar¬ 
row integer quantum Hall strips of opposite chirality if 
ki = 0 and k 2 = ks. This can be motivated by not¬ 
ing that the dominant tunnelings between the inner and 
outer wires induce large gaps for the counterpropagat- 
ing modes at k 2 = fe, and —k 2 = — fe, such that the 
central wire X (Y) forms the bulk of a three wire wide 
integer quantum Hall strip. A right-moving edge mode 
of momentum kx ^ ki = ^ lives in wire A (D), while a 
left-moving edge mode is located in wire B (C). In an 
isolated building block, these edge modes acquire a small 
gap due to both their overlap across the central wire, and 
the direct tunneling between the outer wires. 

With electron-electron interactions, and if the filling is 
reduced, the building blocks can enter fractional quan¬ 
tum Hall statesA Laughlin state at an effective fill¬ 
ing factor V = l/(2m + 1) with positive integer m can 
for instance arise for ki = mk^, and k 2 = {m l)ks. 
While single-particle tunneling between counterpropa- 
gating modes is then forbidden by momentum conserva¬ 
tion, correlated tunnelings can drive the building blocks 
into fractional quantum Hall states. The central wire 
X (Y) then constitutes the gapped bulk of a now frac¬ 
tional quantum Hall strip, while the outer wires again 
host chiral edge modes. A more detailed discussion of the 
physics within a given building block in terms of coupled 
Luttinger liquids is given in Sec. IV Al below. 


To realize 3D phases, I consider the periodic array 
shown in Fig. O I require the tunnelings between the 
inner and outer wires in each building block to be the 
largest energy scale after the bandwidth of the wires, 
and the chemical potential. Additional tunnelings be¬ 
tween wires X and C, and X and D within a unit cell, as 
well as between wires Y and A, and Y and B in neigh¬ 
boring unit cells, which compete with the dominant intra 
building block tunnelings, can then be neglected, and the 
low-energy physics of the array is fully described by the 
integer or fractional quantum Hall edge modes in wires 
A, B, C, and D, as well as the interactions and tunnelings 
between them. 

Focussing first on the integer quantum Hall case ki = 
0, /c 2 = fe, the low-energy dispersions of the edge modes 
are well approximated by Evpkx, where vp denotes the 
Fermi velocity. Along y, I assume neighboring edge 
modes to be coupled by small alternating hoppings tyi 
and ty 2 , whose strengths can be controlled by the inter 
wire distances. I take these hoppings to be shifted by one 
edge mode in the next (x^y) layer. Along z, I connect 
neighboring edge modes by small tunnelings tzi within 
the unit cell, and G 2 between two adjacent unit cells. 
All tunnel couplings are chosen to be positive. 

Labelling the unit cells by an index piny direction and 
q in z direction, the array is described by the low-energy 
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tor of annihilation operators for edge mode electrons 
with momentum kx in wire A, B, C, and D of unit 
cell (p^q). This Hamiltonian is essentially identical to 
the low-energy description of stacked topological insula¬ 
tors analyzed by Burkov and Balents^^ Following their 
calculation, Eq. m is Fourier transformed to momenta 
—Fjay^z < ky^z < k^jkiy^zi where ay {az) is the unit cell 
distance in y (z) direction. The Fourier transform of 
is where k = {kx^ky^kzY" denotes the 3D mo- 
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mentum. Next, it is useful to perform the gauge trans- 
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formations Cj 

^ e^kyayl 2 ^-ikzazl 2 ^ introduce a pseu¬ 

dospin (7 within the A, B, and C, D subspaces, acted on 
by Pauli matrices (Jx,y,zi as well as a pseudospin r be¬ 
tween these two subspaces. The Hamiltonian can then 
be cast into the form H = 'Hk ^k 


T^k —^pkx kJz'kz T T ^?/2) cosi^kytty (Jx 

- {tyl - ty2) Sili{kyay/2) GyTz 

+ {tzi + tz2) cos{kzaz/2) Tx 

+ {tzi - tz2) ^\-a{kzaz/2) Ty . (6) 

After the canonical transformation Gy -F GyTz, Gz -f 
kkz'kzi 'kx 'kxkkxi 'ky 'kyCFx^ f^e diagoualizatiou of the 
r-sector yields H = diag(Hk,+, Hk,-) ^k 


Hk,± =VFkxCrz - {tyl - ty2) ^in{kyay/2)Gy -h M±(k)cr^ , 

( 7 ) 

where the momentum-dependent mass reads 


M±(k)=± yETE+EiEcosEE 

+ {tyl + ty2) COS{kyay/2) . (S) 


IV. PHASE DIAGRAM OF THE 
NON-INTERACTING MODEL 


Despite being described by the same low-energy Hamil¬ 
tonian as stacked topological insulators, I find that the 



FIG. 3: Phase diagram in the non-interacting (left), and 
interacting (right) case with normal insulating (NI), (frac¬ 
tional) quantum anomalous Hall ((F)QAH), Weyl semimetal 
(WS), and single-surface (fractional) quantum anomalous 
Hall (SS(F)QAH) phases. In a lattice model, the filling would 
be 1/2 (1/6) in the non-interacting (interacting) case. 


array of coupled wires can enter a single-surface quantum 
anomalous Hall (SSQAH) phase which is not present in 
Ref. m. This phase is realized for tyi + ty 2 < \tzi — G 2 I 
and tzi < tz 2 , and has a single quantum Hall layer formed 
by the the topmost (for tyi < ty 2 ) or bottommost (for 
tyl > ty 2 ) layer of wires, see below. Interestingly, the 
ratio tyijty 2 thus provides an experimental knob (“quan¬ 
tum Hall switch”) selecting the surface on which the sin¬ 
gle quantum Hall layer, heralded by its gapless edge state, 
appears. 

All other phases of the array are, however, similar to 
Ref. [H. For tyi + ty 2 < \tzi - G 2 I and tzi > tz 2 , the 
system is a normal insulator (NI). If \tzi — G 2 I E tyi + 
ty 2 E tzi + tz 2 ^ the array forms a Weyl semimetal (WS) 
with two gapless Weyl nodes of opposite chirality at k± = 
(0,0,7r/a;2±/C;2o)^, where kzo = arccos(l - [(t^i +^^ 2 )^ - 
{tzi — tz2)^]/‘^tzitz2)/ciz- Surfaces then have a Fermi arc 
between the projection of the Weyl nodes, and show a 
finite Hall conductivity proportional to the distance of 
the projected Weyl nodes. For tyi-\-ty 2 > tzi^tz 2 ^ finally, 
the system is in a quantum anomalous Hall (QAH) phase, 
and exhibits a Hall conductivity of Gxy = jh per (x, y) 
layer of units cells on surfaces with a normal in the (x, y) 
plane. The phase diagram of the array of wires is shown 
in Fig. [3l I have checked that the bulk phase transitions 
as obtained from the effective low-energy model, and the 
presence of a Weyl semimetal phase qualitatively agree 
with a tight-binding calculation that includes all wires 
and tunnelings (also the subleading tunnelings between 
wires X, C, and D, as well as Y, A, and B). 

Besides a blueprint for the constructivist engineering 
of 3D topological states, the coupled-wire construction 
provides particularly simple visualizations of the bulk 
gapped phases based on their hierarchy of couplings. The 
subleading couplings, which compete with the respective 
dominant coupling, are irrelevant, and the system is adia- 
batically connected to an array in which only the leading 
coupling is present. Along the dominant tunnelings, the 
edge modes of the individual building blocks form closed 
loops and open planes of quantum Hall layers, see Fig. 01 
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SS(F)QAH NI (F)QAH 


FIG. 4: Visualization of the gapped phases for tyi < ty 2 . 
Thick solid lines depict the dominant couplings along which 
closed loops and open planes of integer or fractional quantum 
Hall layers form. Open planes are associated with chiral edge 
modes illustrated by thick lines with arrows. 


In the NI phase, all building blocks connect along closed 
loops, resulting in a full gap. The QAH phase, on the 
other hand, consists of open quantum Hall planes alter¬ 
nating with closed quantum Hall loops. Since each open 
plane has a gapless edge mode, the Hall conductivity is 
indeed axy = jh per {x^y) layer of units cells. In the 
SSQAH phase, all building blocks form trivial quantum 
Hall loops, except for the ones on the topmost or bottom¬ 
most layer (depending on the ratio These form a 

single open quantum Hall plane. The WS, finally, occurs 
if the competing tunnelings along y and 2 ; are of similar 
strength. The system then enters a critical phase with 
gapless states of definite momentum (the Weyl nodes), 
which has no simple real-space picture. 


V. ELECTRON-ELECTRON INTERACTIONS 
AND FRACTIONALIZATION 


While the non-interacting array allowed the construc¬ 
tion and analysis of the interesting SSQAH state, and 
gives simple physical pictures of the bulk gapped phases, 
the true power of coupled-wire constructions lies in the 
description of interacting systems. To tackle these, I lin¬ 
earize the spectrum of each wire n = D, X^Y 

within a given unit cell around the Fermi points at 
A^i, 2 , 3 , and decompose the electron operators into 
right {R) and left (L) moving modes as Cn{x) = 
Rn{x)-\-e'^kFLnX where the respective Fermi 
momentum is kprui and r = R, L. For ki = mks and 
k 2 = (m + l)/c 3 with an integer m > 0, momentum con¬ 
servation forbids single-particle tunneling between coun- 
terpropagating modes. The combination of electron- 
electron interactions and tunneling, however, still gen¬ 
erates momentum conserving correlated tunnelings be¬ 
tween neighboring wires. 

In the 2mth order of a perturbation theory in a (con¬ 
tact) interaction V, the tunneling Ia-x between wires A 
and X generates for instance a term proportional to 


tA-x v'^'^ 

X ” + H.c. . (9) 

Connecting different Fermi points, this term trivially con¬ 
serves energy. Momentum, on the other hand, is con¬ 
served for {m + 1) kpRx + '^kpRA = (m -f 1) kpLA + 
'mkFLX^ which is fulfilled for the considered ki = mks 
and k 2 = (m + l)/c 3 . Collecting all leading edge mode 
scatterings at this order in V that conserve energy and 
momentum, I find that the bulk of the A-X-B (C-Y-D) 
building blocks can be fully gapped by 

pt pt m+l 7- mr ^^+1 i rj ^ 

^A{D) -^X(Y) -^X(Y) -^A(D) T rl.C. , (^iUaj 

^X(Y)^^B(C)^~^ ^B(C)^^X(Y)^~^ +H.C. . (10b) 

The leading couplings between the edge modes of the 
building blocks inside a unit cell, on the other hand, are 

+H.C. , (11a) 

+ H.c. , (11b) 

+ H.c. , (11c) 

+ H.c. (lid) 

Couplings between different unit cells will be treated in 
section IV Bl 

A. Luttinger liquid description of an individual 
building blocks 

To understand how these couplings can lead to gapped 
Laughlin states at filling z/ = l/(2m + 1) within each 
building block, it is most convenient to switch to a 
bosonized languageSince the A-X-B and C-Y-D blocks 
can be treated in an analogous fashion, I only detail 
the description of the former. The relevant low-energy 
degrees of freedom are thus the chiral modes close to 
the Fermi points at momentum k^ = ±A^i, 2,3 in wires 
n = A^X^B^ which are bosonized as 

r„(x) = -^e-^^"(-), (12) 

Y Zna 

where is a Klein factor, while a~^ denotes a high 
momentum cutoff. The bosonized fields obey 

[^rn{x),^r'n'{x')\ = 5rr'^nnAxr ^gll{x - x'). (13) 

It is furthermore helpful to also introduce the fields 

+ l)^fln - m^Ln , (14a) 

. (14b) 
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These obey the commutator 

= Srr'Snn'{‘^rn+l)iTrrsgii{x - x') . 

( 15 ) 

Dropping the Klein factors for notational convenience, 
Eq. m translates to sine-Gordon terms 

CO. (SS - SS) , (16a) 

cos (SS - *S) ■ (16b) 

To show that these couplings can stabilize a three-wire 
wide quantum Hall state, I follow Ref. [sH, and note that 
the argument of each sine-Gordon term commutes with 
itself at different positions. These cosines can thus or¬ 
der individually by pinning their arguments to one of 
their minima. In addition, they can order simultane¬ 
ously, since their arguments also commute between each 
other. Because one can always find a Hamiltonian that 
renders the different sine-Gordon terms relevant in the 
renormalization group (RG) sensed it is always possi¬ 
ble to reach this situation at low energies. The result¬ 
ing bulk gapped state has quasiparticle excitations of 
charge e/(2m + 1), as can be inferred from kinks in the 
cosines^i^ and chiral edge modes ^ra and ^lb that 
obey the commutator given in Eq. m, as expected for 
a Laughlin state at filling u = l/(2m + 1)^^ It has fur¬ 
thermore been shown that also the Ghern-Simons term 
associated with a 2D (fractional) quantum Hall effect is 
contained in the coupled-wire construction used for the 
individual building blocks 4^ All of these signatures com¬ 
bined allow to positively identify the gapped phase re¬ 
sulting from the coupling in Eq. (pT|) as an integer (for 
m = 0) or fractional (for integer m > 0) quantum Hall 
state. The integer m, which technically relates to the 
order in perturbation theory in the interaction V that 
generates the required sine-Gordon term, thus physically 
indicates the effective filling fraction u = l/{2m -h 1) of 
the individual building blocks. 

In a larger and isolated building block, where there is 
neither overlap nor direct tunneling between them, the 
edge modes and would be gapless (they simply 
do not appear in the sine-Gordon part of the Hamilto¬ 
nian, and therefore only carry a kinetic energy ^ :^VFkx)- 
In the narrow building blocks considered here, both a fi¬ 
nite overlap through the bulk (the wire X), and a direct 
tunneling between the edge modes exist, such that the 
edge modes are gapped. The leading edge mode tunnel¬ 
ing within the A—X—B building block that is compatible 
with the bulk sine-Gordon terms in Eq. dm is given in 
Eq. ()lla|) . Its bosonized form 

co.(JK-*S) ( 17 ) 

shows that this term can order, and is indeed compati¬ 
ble with the bulk cosines: its argument commutes with 


itself at different positions, and with the arguments of 
the cosines given in Eq. dm. The couplings between the 
edges of the the C — Y — D block and between the two 
building blocks of a unit cell, given in Eq. (ITThli-dTTdli. 
have the bosonized expressions 

cos (iS - iS) . (18a) 

cos (SS - «ic’) . ( 181 .) 

cos (i<S> - ii”>) . (18c) 

The arguments of the individual cosines commute with 
themselves at different positions, and with the arguments 
of the bulk sine-Gordon terms. The arguments of differ¬ 
ent sine-Gordon terms involving the same edge mode do, 
however, not commute, such that these edge mode cou¬ 
plings cannot order simultaneously. 


B. Luttinger liquid description of the 
three-dimensional system 

As in the non-interacting case, I furthermore subject 
the edge modes of the building blocks in the 3D system 
of coupled quantum wires to inter block couplings. I find 
that the leading edge mode couplings compatible with 
the bulk couplings of the individual building blocks in 
Eq. m read 

+ H-c. , ( 19 a) 

+ H.c. , ( 19 b) 

+ H.C. , ( 19 c) 

4 pr^k', + H-c. , ( 19 d) 

where {pq) and {p'q') are the (x^)-indices of neighboring 
unit cells. Like the intra unit cell couplings, these terms 
can be generated by a perturbation theory in an inter¬ 
action E, and inter wire tunneling (the intra unit cell 
couplings in Eq. m can be recovered by setting p = p' 
and q = q'). The couplings in Eq. m conserve energy 
and momentum for the considered case of ki = m ks and 
k 2 = {m ^ 1) ks. Their bosonized form is 


cos 1 

\ARAp'q' 

^LBpq ) > 

(20a) 

cos 1 

ARDp'q' 

- 'i 

^LCpq j > 

(20b) 

cos 

y^RAp'q' 

^LCpqJ ’ 

(20c) 

cos ( 

ARDp'q' 

^LBpq ) > 

(20d) 


where is the bosonized field associated with the 

edge mode r = R, L in wire n = A^B^C^D of unit cell 
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{pq) defined in analogy to Eq. (p!4|) as 

^^Rnpq ~ ^)^Rnpq ~ Lnpq 5 ( 21 a) 

~ ^ ^)^Lnpq ~ Rkpq • (^Ib) 

These fields obey the commutator 

[^^pq{x)^^^^,^,p,q,[x )] =Srr'^nn'^pp'^qq' (2m + 1) 

X i7rrsgn(x — x') . (22) 

The argument of each sine-Gordon term commutes with 
itself at different positions, and can thus order. Like 
within a unit cell, the arguments of different sine-Gordon 
terms involving the same edge mode do not commute, 
and these sine-Gordon terms cannot order simultane¬ 
ously. They do, however, commute with the arguments 
of the sine-Gordon terms stabilizing the integer and frac¬ 
tional quantum Hall states within each building block. 
The couplings in Eq. (EHl) are thus a set of competing 
edge mode couplings that leave the bulk gaps of the in¬ 
dividual building blocks untouched. 

VI. PHASE DIAGRAM OF THE INTERACTION 
MODEL 

Within each A-X-B and C-Y-D building block, where 
all sine-Gordon terms can order simultaneously, a bulk 
gap results from the pinning of the arguments of the 
cosines. For the competing edge mode couplings, how¬ 
ever, the situation is more delicate: not all sine-Gordon 
terms can pin their arguments simultaneously. If the 
couplings have a clear hierarchy in strength, the leading 
sine-Gordon term orders, while the others are suppressed. 
The system is then in a phase that is adiabatically con¬ 
nected to the array being only subject to the leading cou¬ 
pling, while all other couplings vanish. If there is no clear 
hierarchy, the different sine-Gordon terms compete, and 
there is no simple description in terms of pinned fields. 

A. Clear hierarchy of couplings 

Like in the non-interacting case, I take the couplings 
of the inner and outer wires in each building block to be 
the largest correlated tunnelings, which puts each build¬ 
ing block in a fractional quantum Hall state. To analyze 
the effect of the competing edge mode couplings, I study 
an array with the same relative coupling strengths as be¬ 
fore (alternating couplings along y, shifted by one edge 
mode in the next (x, y) layer, and alternating couplings 
along z). If the intra unit cell couplings (113) and (Ed]) 
generalizing the hoppings tzi in the fractional case are 
much stronger than the other edge mode couplings, the 
A-X-B and C-Y-D building blocks pair up within the unit 
cells, and the array is in a normal insulating state. To see 


this, I use that the array is then adiabatically connected 
to a system where the edge modes are only connected by 
this dominant hopping. As shown in the middle panel 
of Fig. m the fractional quantum Hall building blocks 
then form small loops, which indeed corresponds to a 
normal insulating state. If, however, the edge modes 
dominantly pair up along z between neighboring unit 
cells, a single-surface fractional quantum anomalous Hall 
(SSFQAH) phase with a single fractional quantum Hall 
layer on either the top or bottom surface is realized, as 
can be seen in the left panel of Fig. HI If, finally, one of 
the couplings along y is strongest, I obtain a fractional 
quantum anomalous Hall (FQAH) phase analog to the 
above QAH phase, which has one open fractional quan¬ 
tum Hall plane per (x^y) layer of unit cells, see right 
panel of Fig. HJ This phase shows a fractional Hall con¬ 
ductivity of axy = ve^jh per [x^y) layer of units cells, 
where v — ll{2m +1). All of these bulk gapped phases 
are somewhat similar to weak topological insulators in 
that the individual fractional quantum Hall loops and 
planes support gapped quasiparticles of fractional charge 
and statistics confined to their respective 2D subsystem. 
The phase diagram of the interacting array in Fig. [3] lo¬ 
cates the discussed phases in the respective limits of the 
parameter space. 


B. Solution of Luther-Emery type at special points 
in the parameter space 

Most interesting is the fate of the Weyl semimetal. I 
argued that this phase emerges in the non-interacting 
case when several competing tunneling couplings are of 
comparable strength. While the Weyl semimetal phase 
can be identified straightforwardly in a fermionic lan¬ 
guage, as was done in Sec. lEl it corresponds to a com¬ 
plicated regime of competing sine-Gordon terms of com¬ 
parable strength in the Luttinger liquid description, see 
Sec. IV Bl with m = 0. In the interacting array, I find 
that the same competition can also exist between the 
now more complicated sine-Gordon terms. The strik¬ 
ing mathematical analogy of the Weyl semimetal regime 
and the regime of competing edge mode couplings with 
comparable strength in the interacting case hints at the 
possible existence of an extended critical phase separat¬ 
ing the bulk gapped fractional phases. Some support 
for the generic existence of critical phases in 3D comes 
from the superconducting version of the non-interacting 
Hamiltonian in Eq. dlD, which has a gapless supercon¬ 
ducting Weyl phased similar to ^He-A^^ Also spin liq¬ 
uids can exhibit a Weyl phase 

As a first step, it is instructive to analyze the special 
points in parameter space where only one of the compet¬ 
ing sine-Gordon terms is present. For these situations, an 
exact solution of the low-energy Hamiltonian in the spirit 
of a Luther-Emery point is possible for special choices of 
the density-density interactions. For concreteness, let me 
take only the coupling generalizing tz^ to be finite. 
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The Hamiltonian then reads 


= /^quad + <7^ (23) 



Here, i7quad contains the gapless motion along the x di¬ 
rection, and density-density interactions quadratic in the 
bosonized fields. i7quad is thus composed of terms pro¬ 
portional to f dx {dx I ilow intro- 

duce the new fields 


^LCpqi^) 

^RApqi^) 


m -h 1 
2m + 1 

m + 1 
2m T 1 

m -h 1 
2m T 1 

m 1 
2m + 1 


- 

^kDpc2(^) — 


]^^LCpq\^) ’ 

(25a) 

(25b) 

("ri 

(25c) 


These obey the canonical commutator 


[^l^lqix)A^A!’p’g'i^')] =Srr'Snn'Spp'Sqq> iwrsgnix - x') . 

(26) 

Rewriting the Hamiltonian in terms of these new degrees 
of freedom yields 



where H^q^ad is still composed of products of linear deriva¬ 
tives of the new bosonized fields ^r^lq. Because the fields 
^r^pq obey the canonical commutator of the bosonized 
fields describing chiral fermions, I can now rewrite the 
Hamiltonian as 


H = + ffint, (29) 

= [dxJ2¥;^'>\x)n^^\x)¥;^\x) , (so) 

^ P,q 


n'r\x) 


/- 

-iudx 

: 0 

Am) 

f^Zi 

0 \ 


0 

iudx 

0 

Am) 

^^Zl 


Am) 

0 

iudx 

0 

v 

0 

Am) 

l^Zi 

0 

-iudx) 


(31) 


^pq 


1 

\/2'Ka 



(m) 

RApq 


e 


— 


(m) 

LBpq g 


—i<^> 


(m) 

LCpq g 


^RDpqy 


(32) 


where i^int describes only the density-density interactions 
between the different modes, and where u is the velocity 
associated with the fields ^^Iq (which is for simplicity 
assumed to be identical for all modes). I can now invert 
the bosonization procedure for the components of 
for all m using the refermionization procedure^^iS Diag¬ 
onalizing the refermionized version of as in Sec. IHII 
then allows to calculate its spectrum for all m. As follows 
from Sec. nni upon setting all tunnel couplings except tz^ 
to zero, the system is in a normal insulating state for 
y 0, and has a critical point for = 0. 

It is interesting to also determine the charge associated 
with the operator , which can be inferred from the 

bosonized expression of the total charge density operator. 


Plot — ^ ^ ^ ^ Rnpq ^Lnpq^ 

pq n=A,B,C,D 

= Y. Y. 27 r( 2 m + 1) 

pq n=A,B,C,D ^ ’ 

= Y Y 2n{2m + 1) ’ 

pq n=A,B,C,D ^ ' 

( 33 ) 

where e denotes the charge of an electron. Commut¬ 
ing ptot and I find that these latter exponen¬ 

tials create quasiparticles of charge e/(2m + 1). This 
finding agrees with the observation that 27r-kinks (the 
smallest possible kinks) in any of the bulk sine-Gordon 
terms given in Eq. (|20|) carry charge e/(2m + 1)^ and 
that quasiparticle operators add local excitations above 
the ground state. Indeed, kinks in the bulk sine-Gordon 
terms precisely correspond to local excitations above the 
ground state. Since the analysis performed here for only 
being finite can be directly transposed to any other 
coupling, I find that if only one of the competing tun¬ 
nel couplings is present, the system can be described in 
terms of interacting fermionic quasiparticles of fractional 
charge e/(2m + 1) that have a critical point, at which 
they disperse linearly in x direction. In the spirit of a 
Luther-Emery point, one can now consider fine-tuning 
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the density-density interactions in the system to reach 
^int = 0. At this point, the non-interacting spectrum 
of is the true spectrum of the system, and the 

strongly interacting electron system has been mapped 
exactly to free fermions of fractional charge. As for 
a one-dimensional Mott insulator at half-filling with a 
Luttinger liquid parameter K = 1/2, these emerging 
fermions are complicated objects that are not to be con¬ 
fused with the original electrons in the system^^ 

C. Competing couplings of comparable strength 

When more than one of the competing sine-Gordon 
couplings is present, the transformation in Eq. (|25|) is 
not sufficient to rewrite the Hamiltonian in terms of (free) 
fermions anymore. For the non-interacting case m = 0, 
I showed in Sec. that the effect of the additional cou¬ 
plings is to turn the critical point associated with tzi = 0 
considered in the last subsection into an extended criti¬ 
cal phase, the Weyl semimetal. Put in simple terms, the 
additional couplings allow the gapless quasiparticles to 
move in 3D. For finite m, one may suspect that the crit¬ 
ical point associated with = 0 considered in the last 
subsection also turns into an extended critical 3D phase 
when the additional couplings are turned on. Similar to 
a discussion supporting fractional topological insulators 
in Ref. [HU, this scenario can be supported by rewriting 
the Hamiltonian as 


where is the number of unit cells in y {z) direc¬ 

tion. 

Within the region of the phase diagram hosting the 
Weyl semimetal in the non-interacting case m = 0, also 
the interacting array has excitations associated 

with momentum ky = 0 and = Tr/a;^ ± kzo that do 
not appear in the Hamiltonian i^cos^ but only in ii^quad 
(which accounts for the possibly interaction-renormalized 
kinetic energy along the x direction). For m = 0, these 
are precisely the gapless electronic quasiparticles at the 
Weyl nodes. Outside the parameter regime hosting the 
Weyl semimetal in the non-interacting case, the Hamilto¬ 
nian matrix 'Hkykzi^) does not vanishes for any {ky^kz). 
As a consequence, all excitations create kinks of non¬ 
trivial energy in the sine-Gordon Hamiltonian i^cos, and 
are thus gapped. This suggests that the topology of the 
phase diagram is the same in the interacting and non¬ 
interacting cases. 

As discussed in the last subsection, the fermionic char¬ 
acter of the operators for m = 0 can most easily be 
derived by refermionizing its components For 

m > 0, however, the fields do not obey the same 

commutator as the non-interacting fields ^rnpq. This im¬ 
plies that the refermionization procedure does not carry 
over to the interacting case, and that the simple opera¬ 
tors do, for general m, not satisfy the canonical 

fermionic commutation relations. This is in agreement 
with the observation that the operators 


H — -Hquad T 
Hcos ^ J 


cos 7 


p,q,p'q' 


ypp'qq' _ I ^11 i^) Hu 


(34) 

, (35) 

pp'qq\x)\ 


(Hll «« {x) H 


(36) 


22 


\I/(”») = __ 


(e' 


^RApq ^ g ^^LBpg ^ g ^ g ^^RDpq )"^ 

(37) 


Here, the (2x2) matrices (x) contain the pref¬ 

actors of the sine-Gordon terms. By adjusting the tun¬ 
nel couplings and interaction strengths, can be 

tuned to be of the same form for all integers m, in¬ 
cluding the non-interacting case m = 0. Gonsequently, 
the same transformations that diagonalize the Hamilto¬ 
nian matrix for m = 0 also do so for positive 

integers m. These transformations follow from Eq. m 
and its discussion in Sec. upon setting vpkx 0, 
and replacing the non-interacting tunnelings ty^z, 1,2 by 
the prefactors of the corresponding sine-Gordon terms. 
As one important step in this diagonalization, I high¬ 
light the Fourier transformation along y and 2 ), which 

yields Fcos = 
with T 


(t) = 1 




f i$(0) f , A™ /oo^ 

Q^^rnpq = I | I g {-r)npq | ( 3 o) 

correspond to products of 2m + 1 fermionic operators, 
rather than a single canonical fermion. The operators 
(x) with m > 0 thus do not correspond to indi¬ 
vidual 3D quasiparticle excitations of the array of wires. 
Further extensions of the coupled wire construction are 
needed to determine if also the 3D quasiparticles at fi¬ 
nite m have linearly dispersing band touchings similar 
to Weyl nodes. Similarly, the properties of the interact¬ 
ing analog of the Weyl semimetal phase, such as its re¬ 
sponse to electro-magnetic fields, cannot be resolved by 
the present analysis based on a language of competing 
sine-Gordon terms. 

Despite the fact that the present bosonized calculation 
does not allow to access the true quasiparticles in the 
generalized Weyl semimetal regime for finite m, I judge 
the striking analogy between the interacting and non¬ 
interacting cases, and (for all m) the existence of oper¬ 
ators that do not appear in i^cos in tbe regime hosting 
the Weyl semimetal for m = 0, and only in this regime, 
to strongly support the existence of an extended criti¬ 
cal phase separating the bulk gapped fractional phases 
emerging from the critical points discussed in Sec. IVIBi 
similar to the Weyl semimetal emerging from the analo¬ 
gous points point for m = 0. While one option for a pos¬ 
sible critical phase would be a standard Weyl semimetal. 
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I note that transitions between the bulk gapped phases 
correspond to a rearrangement of only the edge states of 
the individual fractional quantum Hall building blocks. 
Their bulk gaps in the X and Y wires, however, never 
close. I thus speculate that the quantum phase transi¬ 
tions of the interacting array, and the possible gapless 
critical phase inherit some of the exotic properties of 
fractional quantum Hall edge modes. A critical phase 
could for instance be composed of fractionally charged 
fermionic 3D quasiparticles (somewhat related ideas are 
discussed in Refs. This scenario is supported by 

the quasiparticles at the Luther-Emery-type points con¬ 
sidered in Sec. IVIBl 


VII. SUMMARY AND CONCLUSIONS 

In this work, the coupled-wire approach has been 
adapted for the analysis of interacting topological phases 
in 3D. More precisely, I analyzed the 3D system resulting 
from connecting narrow 2D integer and fractional quan¬ 
tum Hall building blocks along different directions. The 
coupled-wire language is thereby used for the microscopic 
description of the individual building blocks, see Sec. HIl 
and the inter block couplings. For the non-interacting 
array of wires detailed in Sec. imi I found that the phase 
diagram includes a normal insulating, a quantum anoma¬ 
lous Hall, a Weyl semimetal, and a single-surface quan¬ 
tum anomalous Hall phase, see Sec. IIVI The surface on 
which this latter phase appears can be controlled by the 
ratio of two tunnel couplings. The coupled-wire language 
offers simple real space visualizations of the bulk gapped 
phases in terms of closed loops and open planes of 2D 
quantum Hall subsystems. With electron-electron inter¬ 
actions, I found that the array can be in analogous bulk 
gapped phases composed of closed loops and open planes 
of fractional quantum Hall states, see Sec. [V] and Sec. lVIl 
Some of these phases have topologically protected frac¬ 


tional edge states. In analogy to the non-interacting Weyl 
semimetal, the bulk gapped phases are separated by a 
regime in which certain excitation operators do not ap¬ 
pear in the sine-Gordon part of the Hamiltonian. This 
supports the existence of an exotic critical phase simi¬ 
lar to a Weyl phase. The phase transitions of the ar¬ 
ray are of exotic nature, and correspond to a rearrange¬ 
ment of fractional quantum Hall edge modes. If only 
a single sine-Gordon term is present, the exotic nature 
of the phases and phase transitions could be revealed 
by an exact mapping, which identified the quasiparticles 
at special Luther-Emery points as fermions of fractional 
charge e/(2m +1). The properties of the critical regime 
outside these special points will be addressed in future 
work. Other interesting directions include the use of chi¬ 
ral spin liquids, topological superconductors, and states 
with a more complex edge structure (such as Moore-Read 
states) as 2D building blocks, as well as the coupled-wire 
analysis of non-Abelian 3D phases with string-like ex¬ 
citations. Experimentally, the proposed array could be 
realized with cold atoms, in which the necessary ingre¬ 
dients (fermions subject to spin-orbit coupling and mag¬ 
netic field^i^ complex lattices^ and superlattices host¬ 
ing for instance even resonating valence-bond states^) 
have been demonstrated. 

Note added. Recetly, a related preprint discussing 
coupled-wire constructions of 3D integer and fractional 
topological insulators appeared^^ 
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